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C$ ■ Abstract. We study the difference between thermally produced fermionic and 

bosonic hot dark matter in detail. In the linear regime of structure formation, their 
distinct free-streaming behaviours can lead to pronounced differences in the matter 
power spectrum. While not detectable with current cosmological data, such differences 
will be clearly observable with upcoming large scale weak lcnsing surveys for particles 
as light as whdm ~ 0.2 eV. In the nonlinear regime, bosonic hot dark matter is not 
subject to the same phase space constraints that severely limit the amount of fermionic 
hot dark matter infall into cold dark matter halos. Consequently, the overdensities in 
fermionic and bosonic hot dark matter of equal particle mass can differ by more than 
a factor of five in the central part of a halo. However, this unique manifestation 
of quantum statistics may prove very difficult to detect unless the mass of the hot 
dark matter particle and its decoupling temperature fall within a very narrow window, 
1 < "T-HDivi/eV < 4 and g» < 30. In this case, hot dark matter infall may have some 
observable consequences for the nonlinear power spectrum and hence the weak lensing 
convergence power spectrum at I ~ 10 3 — > 10 4 at the percent level. 
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1. Introduction 

With the advent of precision measurements of the cosmic microwave background (CMB), 
large scale structure (LSS) of galaxies, and distant type la supernovae, a new paradigm 
of cosmology has been established. In this new standard model, the geometry is flat 
so that Qtotai — 1) an d the total energy density is made up of matter (Q m ~ 0.3) 
[comprising of baryons (Qb ~ 0.05) and cold dark matter (^cdm ~ 0.25)], and dark 
energy (fix ~ 0.7). With only a few free parameters this model provides an excellent 
fit to all current observations PjJ 121 ED E] ■ In turn, this allows for constraints on other, 
nonstandard cosmological parameters. One very interesting possibility discussed widely 
in the literature is a subdominant contribution to the total energy density in the form 
of neutrino hot dark matter (HDM) and hence limits on the neutrino mass. 

In general, any fermionic dark matter species if which decouples while relativistic 
obeys the relation 

„ , 2 10-75 

g*^ 183 ev 

where is the %fj particle's mass, g its number of internal degrees of freedom, and g* t ^ 
the effective number of degrees of freedom in the plasma at the time of ip decoupling. 
By entropy conservation, g^ can be related to the effective temperature of the species 
via g^^T^ = g* )U T^, with the subscript v referring to the standard model neutrinos. 
Standard model neutrinos decouple in the early universe at a temperature of order 
2^3 MeV, when g* = 10.75. Thus, just from demanding Q v < fi total — 1, one finds 
the well known upper limit on the neutrino mass [SHU], vn v < 46 /N v eV, assuming N v 
neutrino flavours with degenerate masses. 

However, a much stronger limit on the neutrino mass can be derived by noticing 
that the thermal history of HDM is very different from that of CDM. By definition 
HDM becomes nonrelativistic only at very late times. At early times free-streaming 
of the HDM particles causes essentially all of their own perturbations to be erased on 
scales below the free-streaming length [Jj, 
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leaving only perturbations in the nonrelativistic matter (CDM and baryons). On scales 
larger than Aps, however, HDM behaves like CDM. Thus, the net result is a suppression 
of the overall level of fluctuations on scales below Afs- 

In terms of the present matter power spectrum, P(k,r) = \5\ 2 (k,r), where S(k,r) 
is the Fourier transform of the density perturbations S(x,t), early free-streaming leads 
to a suppression of power at k ^ 27r/Aps by roughly 

P{k,T) il m 

where the factor eight is derived from a numerical solution of the Boltzmann equation 
jUJ- Equation (jl.3j) applies only when Q m ; when dominates, the spectrum 
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suppression becomes exponential as in the pure HDM model. Note that the free- 
streaming scale becomes very small for very large values of g* (~ 10 3 ). In this case, the 
species should be better known as warm dark matter (see [B] for a recent discussion of 
mass bounds on warm dark matter). 

Of course, a suppression of small scale power in the present day matter power 
spectrum could have its roots in a number of other factors, such as a lower matter density, 
a higher relativistic energy density, a primordial power spectrum with broken scale 
invariance, etc.. Therefore, in order to constrain the HDM energy density and hence 
the mass of the HDM particle, it is necessary to combine LSS with CMB measurements. 
Currently available cosmo logical data constrain the sum of light, standard model 
neutrino masses to J2 m u — 0.5 — ► 2 eV, depending both on the data sets used and 
on assumptions about other cosmological parameters [HI EH HU Il2l HH1 UH ESI EH EZj 

While cosmological limits on HDM have mainly been discussed for neutrinos, they 
apply equally well to bosons and thermal relics with g* > 10.75 jTHj. Many such 
hypothetical HDM candidates exist, including the likes of axions jTHl ED], majorons, 
gravitinos [21], axinos [22J, etc.. An even more exotic scenario recently discussed in the 
literature is the possibility that neutrinos violate Fermi-Dirac spin statistics and behave 
like bosons [231 121] ■ It is clearly worthwhile to study the differences between these 
various HDM species arising from their quantum statistics. As we shall see later, there 
are fundamental differences in their clustering properties which might be observable 
with future, high precision probes. 

Indeed, the purpose of this paper is to study the differences between fermionic and 
bosonic HDM. From here on we consider two generic HDM species: A single massive 
Majorana fermion ip with g = 2, and a single massive scalar <fr with g — 1. Unless 
otherwise indicated, we shall assume, throughout the present work, a cosmological 
constant A, three massless standard model neutrinos, and the following cosmological 
parameters: {fl m , fl\, h, as} = {0.3, 0.7, 0.05, 0.7, 0.9}. 

The paper is structured as follows. In the next section we discuss structure 
formation in the linear regime. Section El deals with nonlinear structure formation, 
including hot dark matter infall into CDM halos. We consider gravitational lensing as 
a probe of structure formation in section HJ and section El contains our conclusions. 



2. Structure formation in the linear regime 

A thermal relic that decouples from the cosmic plasma while relativistic assumes either 
the relativistic Fermi-Dirac (FD, "+") or Bose-Einstein (BE, "— ") distribution, 

expyo/ i J ± 1 

If the particle is stable and no heavier particle decays into it, its phase space distribution 
is preserved with the expansion of the universe. Today, a single massive g = 2 Majorana 
fermion if) and a g = 1 scalar <fi are expected to contribute to the total energy density 
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Figure 1. Linear power spectra for two different AHCDM models. The blue (dotted) 
line shows a model with three massless neutrinos and one massive Majorana fermion, 
contributing fl^ — 0.02. The red (solid) line shows the same, but with a massive scalar 
instead. The black (dashed) line is the standard ACDM model with no HDM. Note 
that these spectra have been normalised to have the same amplitude on large scales. 



at fractions of 



Q = 2 x — x (2.2) 
* 183 eV 1 ; 

¥^lx^x^L. (2 .3) 
3 183 ev 

Already, we can see that for any given f^DM; the mass of the scalar is necessarily a 
factor 3/2 higher than that of the fermion. 

With regard to structure formation, the amount of suppression in the matter power 
spectrum due to HDM is essentially dependent only on the ratio f2 H DM/^m, i-e., it 
depends only on the lack of dark matter. Therefore, on both very small and very large 
scales, fermion and scalar models with identical Qhdm are virtually indistinguishable 
to the naked eye. However, for scales close to the free-streaming length (and therefore 
also to the scale of matter-radiation equality for sub-eV hot dark matter), there is a 
pronounced difference. The reason is that for the same Ohdmj the more massive scalar 
particle exhibits less free-streaming. This in turn allows for an increased fluctuation 
amplitude around the free-streaming scale. In Figure we show the power spectra for 
two models with identical Qhdm = 0.02 and g* = 10.75 (corresponding to = 0.9 
eV for the fermion model and = 1.3 eV for the scalar model). As can be seen, the 
scalar model has more power at k ~ 0.01 — ► 0.1 h Mpc. For comparison we also show 
the standard ACDM model, which has the same parameters, except that Qhdm = 0. 
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Figure 2. Number densities of HDM expected for — 10.75 and the indicated halo 
and HDM masses at z = 0, normalised to the cosmological average hi. Red (solid) 
lines represent bosons, while blue (dotted) lines denote fermions. These densities are 
obtained from solving the Vlasov equation using the method of assuming CDM 
halos with the Navarro-Frenk- White density profile [equation l|B.8(l ], 

3. Nonlinear structure formation 

At late times (z < 100), fluctuations in the mass density field can become larger 
than unity. Linear perturbation theory breaks down, and structure formation enters 
a nonlinear phase with the collapse of overdense regions into gravitationally bound 
objects. In hierarchical CDM cosmologies, nonlinear gravitational collapse begins at 
small scales, and the nature of the CDM particle determines the size of the first objects. 
Subsequent mergers of these small systems give rise to larger structures — from the halos 
of dwarf galaxies to the filaments and the voids. In terms of Fourier decomposition, the 
transition from the linear to the nonlinear regime in the density field occurs roughly at 
k ~ 0.2 h Mpc -1 today. 

The effect of HDM on nonlinear clustering and its statistics is twofold. Firstly, 
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Figure 3. Left: Nonlinear power spectra at z — for g» = 10.75 and various HDM 
masses, relative to the standard ACDM power spectrum. Solid lines denote bosons 
of masses (top to bottom) {0.3,0.6,1.2} eV, dotted lines are for {0.2,0.4,0.8} eV 
fermions, while dashed lines represent {0.3,0.6,1.2} eV fermions. These masses are 
chosen such that they satisfy either fi^, = fi^ or iriw, = m^. Note that all spectra here 
have been normalised to a fixed value of erg. Right: Relative contributions from HDM 
infall for the same HDM scenarios. 



free-streaming of HDM in the early (linear) stages of structure formation suppresses 
density fluctuations on scales below Afs- This directly limits the number of small scale 
structures that can be formed subsequently in the nonlinear phase. 

Secondly, at late times, secondary infall of HDM into existing CDM halos becomes 
possible when the former's mean velocity, 

V 9* J Vm H DM/ { 6 BE 

drops below the velocity dispersion of the astrophysical system concerned. A typical 
galaxy cluster (M ~ 10 14 — > 10 15 M©) has a velocity dispersion of about 1000 km s _1 
today; a typical galaxy (M ~ 10 12 M ), about 200 km s~ . Thus, for whom ~ 1 eV, 
a good fraction of HDM particles can be expected to reside presently in halos across a 
wide mass range. Figure El shows a sample of HDM overdensities expected for a scalar 
and a Majorana fermion for a range of halo masses. These are obtained from solving 
the Vlasov equation (jA.l|) . assuming CDM halo density profiles of the Navarro-Frenk- 
White form [equation (|B.8|) ] (2E1EE!- A brief description of this calculation can be found 
Appendix ~£| For more details, see reference |27j . 



in 



Observe in Figure El that at the sub- Mpc scales, bosons cluster considerably more 
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Figure 4. Nonlinear power spectra at z = for two different values of g*, relative to 
the standard ACDM power spectrum. Red (dark) lines denote fermions, while green 
(light) lines represent bosons. Solid lines show the power spectra with HDM infall, 
dashed lines show them without. 

efficiently than do fermions. This is because, at any given temperature, the unperturbed 
phase space distribution for bosons is more skewed towards the low momenta, making 
the particles more susceptible to capture by the CDM halos. Furthermore, the final 
HDM phase space distribution in a halo is, in general, subject to constraints imposed by 
Liouville's theorem, and hence must not exceed the maximum of the initial unperturbed 
distribution function, / final < /£jtf a i EE]- For fermions with no chemical potential, the 
maximum /Stfai = 1/2 occurs at p — 0. Using this line of argument, Tremaine and 
Gunn have famously derived a bound on the maximum possible neutrino density in a 
halo [2H1 ED!- Bosons, however, are technically not subject to the same constraint, since 
no finite maximum exists in their initial phase space distribution. Thus, correspondingly 
in Figure the overdensity of bosons is seen to be much higher than that of fermions 
with the same mass. 

Note, however, that the Tremaine-Gunn bound is still applicable to bosons in a 
statistical sense because only a very small fraction of all thermally produced bosons 
reside in low momentum states [HH E2|- (One notable exception is axion cold dark 
matter which consists of low mass bosons in a condensate.) 

3.1. The nonlinear power spectrum 

In Figures El to we show predictions for the nonlinear power spectrum at z = for a 
variety of scenarios involving bosonic and fermionic HDM. These spectra are calculated 
from a phenomenological approach known as the halo model EU ESI EE] , extended 
to include contributions from HDM infall jSZj- While the halo model will not serve for 
high precision cosmology, it suffices to illustrate the typical scales of the various effects 



Is it possible to tell the difference between fermionic and bosonic hot dark matter? 8 




Q5 I i i i i i i i i i i i i i 1 1 1 1 i i i i i 1 1 1 

10" 1 1 10 10 2 10 3 

k [h Mpc" 1 ] 



Figure 5. Nonlinear power spectra at z — for fermions with various values of 
and but the same Q^h 2 — 0.011, relative to the standard ACDM power 
spectrum. Solid lines show the power spectra with HDM infall, dashed lines 
without. From top to bottom: — {1,1.41,2.00,2.82} eV, corresponding to 

g* = {10.75,15.16,21.5,30.32}. 



due to HDM. See Appendix B for more details. 

As is evident in the Figures, early HDM free-streaming is generally the dominant 
influence on the shape of the nonlinear power spectrum. Ignoring contributions from 
late-time infall, one can see that early free-streaming alone already distinguishes between 
fermions and bosons, and between cases with different g*'s at the percent level. Late- 
time HDM infall causes additional losses of power at k > 1 h Mpc -1 , because HDM 
clustering tends to "smooth out" the overall matter distribution within the individual 
halo. For the = 10.75 cases (Figure EJ), losses from infall are typically no more 
than one percent, and occur at scales corresponding to roughly ten times the current 
free-streaming wavenumber of the particle kps, 



k FS - / x 0.11 (_M 1/8 (^L) h Mpc -1 . (3.2) 



1 + z V 10.75/ V eV / 
The analysis of [37] also finds a reduction of power due to infall at the 0.1 — > 1 % level, 
albeit at a scale of k ~ 0.5 h Mpc -1 . This discrepancy could be due to the linear method 
used in their clustering calculations. 

Infall increases in importance for cases with g* > 10.75. This is because these 
particles are colder, and for the same Qhdm, are heavier; both attributes allow them to 
cluster more efficiently Indeed, as shown in Figures 0] and HDM infall in any one 
scenario can alter the resulting power spectrum to the extent that it mimics another 
scenario, and is therefore potentially important for high precision cosmology. Examples 
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of thermal relics with large g*'s include very light gravitinos (g* ~ 100, corresponding to 
a decoupling temperature of a few GeV to a few TeV), and thermally produced axions 
(10.75 < g* < 80, corresponding to decoupling during the QCD phase transition). 

4. Observational probes 

4-1. Weak gravitational lensing 

Weak gravitational lensing of distant galaxies by intervening large scale structure 
provides a unique method to directly map the matter distribution in the universe 
(e.g., |HE])- Perturbations in the gravitational potential along the line of sight induce 
distortions in the images of distant galaxies. The images may be magnified (convergence) 
and/or stretched (shear). For weak lensing of a large ensemble of sources and lenses, it 
can be shown that the convergence and the shear have identical statistical properties. 
Henceforth, we shall consider the power spectrum of the convergence as representative 
of the lensing features. 

The convergence is essentially an integral over all the deflectors between us and 
the source galaxies, weighted by the source galaxy distribution. Its power spectrum at 
multipole £ can be related to the matter power spectrum P(k, z) via |3Hl EH3 E0| 



in a flat universe. Here, x — Jo dz/H(z) is the comoving radial distance, and Xh is the 
distance to the horizon. The weak lensing weighting function, 



encapsulates information about the source galaxy redshift distribution n(x)- For our 
purposes, n(x) ma Y be taken to be n(z) oc z a exp[— (z/zof], where the parameters zo, a 
and (3 are estimated from existing deep redshift surveys for the weak lensing survey at 
hand jlTJ|321- The prefactor is determined by the normalisation condition / dz n(z) = 1. 
In the following, we shall take z = 1, a = 2 and (3 = 2, typical of most proposed large 
scale lensing surveys |42j . With these parameters, the weak lensing weighting function 
flPD peaks at z ~ 0.4 -> 0.5. 

In any lensing survey, the ability to measure Ce is constrained on large scales by 
a sample variance due to finite sky coverage, and on small scales by the finite number 
of available galaxies. Neglecting non-Gaussian corrections, the statistical error in Ce is 
estimated to be jS3 ED] 



where / s k y is the fraction of the sky covered by the survey, n ga i is the surface density 
of galaxies, and 7 rms is the rms shear per galaxy (from noise and intrinsic ellipticity). 
These parameters are generally survey dependent. Proposed wide-field surveys such as 




(4.1) 




(4.2) 




(4.3) 
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Figure 6. Top: Weak lensing convergence power spectra for the models shown in 
the left panel of Figure The red (dark solid) line denotes fermions, the green (light 
solid) line represents bosons, while the black (dotted) line corresponds to the standard 
ACDM model. Bottom: The same two convergence spectra normalised to the standard 
ACDM convergence spectrum. Solid lines show the spectra with HDM infall, dashed 
lines show them without. The grey shaded area represents errors expected for a LSST- 
type survey, smeared over bands of width £/4. Figured left shows the same plot at 
i = 10 2 — > 10 5 for greater clarity. 



the Supernova/Acceleration Probe (SNAP) |Hj are expected to cover some 1000 deg 2 
(/sky ~ 0.03) of the sky; the coverage of the more ambitious Large Synoptic Survey 
Telescope (LSST) will likely be another tenfold. In the ensuing analysis, we shall 
adopt 7 rms = 0.15, n gal = 50 arcmin -2 , and / s k y = 0.5, parameter values corresponding 
to an optimal LSST-type survey. After smearing over bands of width £/4, we have 
approximately 

jQ^^HLi] (4.4) 
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Figure 7. Weak lensing convergence power spectra for the models shown in Figure 
relative to the standard ACDM convergence spectrum. Red (dark) lines denote 
fermions, while green (light) lines represent bosons. Solid lines show the spectra with 
HDM infall, dashed lines show them without. 

as the effective error |43| . 

Figures El to El show the expected weak lensing convergence power spectra for the 
various HDM scenarios considered in the previous section. The general shapes of the 
convergence spectra are similar to those of their corresponding matter power spectra 
shifted according to k ~ £/1000, and reflect predominantly the early free-streaming 
behaviours of the different HDM species. Contributions from infall, however, are 
considerably reduced (by a factor of 2 — > 3) relative to the z = matter power spectrum. 
This is because weak lensing is primarily sensitive to the z ~ 0.4 — > 0.5 universe, where 
the amount of HDM infall is correspondingly lower. 

In Figures 13 and |Hl we see that the contributions from infall are typically at the sub- 
percent level at I ~ 10 3 — > 10 4 . As with the nonlinear power spectrum, the contributions 
increase with and mjiDM- However, because the scale at which HDM infall is manifest 
depends on the current free-streaming scale of the HDM particle [cf. equation (J3.2)) ]. the 
signatures of infall are also shifted accordingly to larger values of £, where the statistical 
errors are larger. The net result is that there may only be a very narrow window in tbhdm 
and g* in which HDM infall has some observable consequences in weak lensing, even for 
an optimal LSST-type lensing survey. Judging from Figures and |Hl we estimate this 
window to be 1 < mnDM/eV < 4 and g* < 30. 

Nevertheless, fermionic and bosonic HDM can still be distinguished by their early 
free-streaming behaviours. In Figure El we see that this is possible for an LSST-type 
lensing survey for particles as light as m H DM ~ 0.2 eV. 

4-2. Strong gravitational lensing 

Looking at Figure El it is clear that the main difference between fermionic and bosonic 
HDM is in the central regions of halos. The matter distribution in the central parts 
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Figure 8. Left: Weak lensing convergence power spectra for the £l^h 2 = 0.011 
models shown in Figure relative to the standard ACDM convergence spectrum. 
Solid lines show the spectra with HDM infall, dashed lines without. The masses 
of the HDM particles are (top to bottom) {1,1.41,2.00,2.82} eV, corresponding to 
<7* = {10.75,15.16,21.5,30.32}. Right: Same as the left panel, but zoomed in on 
I = 10 2 -> 10 5 . 



of dark matter halos can in principle be probed by strong gravitational lensing, either 
by very precise measurements of a single lensing system ^01 113 HE] , or by measuring 
a very large sample of systems (for instance using a supernova survey |49|). However, 
even though difference between fermions and bosons is large, the total contribution of 
HDM in the central part of a halo is only a minute fraction of the total matter density. 
The ratio Phdm/pcdm decreases with decreasing r and therefore HDM is most easily 
probed on scales comparable to, or larger than, the virial radius where weak lensing is 
most efficient. 

5. Conclusions 

We have studied the difference between fermionic and bosonic hot dark matter in detail. 
In the linear regime of structure formation, the most important difference is that for 
equal contributions to the present energy density, the mass of a scalar species is higher 
than that of a fermion by a factor 3/2. This in turn means that scalars have a smaller 
free-streaming length and therefore less suppression of fluctuation power. This effect 
is most pronounced around the free-streaming scale, where the difference can be of 
order 30 — > 40 % in terms of the present day linear power spectrum. At present, the 
difference is too small to be discerned by available cosmological data (as also discussed 
in the context of nonthermally distributed neutrinos |SD1), but may be visible in the 
next generation of probes for large scale structure and in proposed wide-field weak 
gravitational lensing surveys. 

In the nonlinear regime, there is a fundamental difference between fermions and 
bosons arising from the fact that the phase space distribution for fermions is subject to 
constraints imposed by Liouville's theorem, while for bosons there is no formal upper 
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Figure 9. Left: Weak lensing convergence power spectra for the = 10.75 models 
shown in Figure |3J relative to the standard ACDM convergence spectrum. Solid 
lines denote bosons of masses (top to bottom) {0.3,0.6, 1.2} eV, dotted lines are for 
{0.2, 0.4, 0.8} eV fermions, while dashed lines correspond to {0.3, 0.6, 1.2} eV fermions. 
Right: Same as the left panel, but zoomed in on I = 10 2 — » 10 5 . 



bound. This difference has a direct impact on the particles' clustering properties. In 
order to probe this phenomenon, we have studied late-time hot dark matter infall into 
cold dark matter halos, and find that, in the central part of a halo, the difference 
between the fermionic and bosonic overdensities can be more than a factor five. While 
this difference appears large, it is not likely to be observable by direct means (e.g., by 
strong gravitational lensing) simply because the density of hot dark matter is so small 
in the central regions of halos. 

Instead, we have considered the implications of hot dark matter for the nonlinear 
power spectrum, k > 0.2 h Mpc -1 . These are the scales at which future weak lensing 
surveys have the most statistical power. Using a halo model approach, we find that 
the early free-streaming behaviour of the hot dark matter species remains the dominant 
influence on the shape of the resulting nonlinear power spectrum. In some cases, notably 
thermal relics with > 10.75, hot dark matter infall can contribute an additional loss of 
power at the percent level at z = 0, depending on their masses and quantum statistics. 
This suggests that infall can potentially play a decisive role in high precision cosmology. 
However, because both the amount of hot dark matter infall and the scale at which it 
is manifest depend crucially on the values of g* and mnDM, signatures of infall may end 
up unobservable in the weak lensing convergence spectrum, even with a Large Synoptic 
Survey Telescope-type lensing survey, unless g* < 30 and 1 < mnDM/eV < 4. 



Is it possible to tell the difference between fermionic and bosonic hot dark matter? 14 

The bottom line is that the difference in the hot dark matter contribution arising 
from quantum statistics is a very interesting example of statistical mechanics, but that 
it may prove to be very difficult to observe. On the other hand, imprints on the matter 
power spectrum from different early free-streaming behaviours extend to both the linear 
and the nonlinear scales. These could well be observable with future experiments if the 
mass of the hot dark matter species is sufficiently large (say, m H DM ^ 0.2 eV). 



Appendix A. Hot dark matter clustering 

The problem of late time HDM clustering in the presence of a CDM halo can be treated 
quantitatively with the nonrelativistic Vlasov equation (e.g., 

^ + — • ^ - a mi V(f) — = 0, i = CDM, HDM, (A.l) 
or anii ox op 

which tracks the evolution of the phase space densities fi(x,p,r) under the influence 
of the peculiar gravitational potential 4>(x, r) as functions of comoving position x, 
momentum p, and conformal time r. In turn, <p(x, r) is related to the local density 
fluctuations 5i(x,r) via the Poisson equation, 

W = — EA*(*,r), (A.2) 
a i 

where 

5i(x,r) = Pl ^ X J _ i pi(x,T) = nii f d 3 p fi(x,p,T), (A.3) 
Pi J 

and pi denotes the ith comoving mean density. 

A number of methods are available for solving equations (jA.ljl and (|A.2j) — from 

numerical simulations to semi-analytical approaches under a variety of approximation 

schemes — with varying degrees of accuracy. In this work, we use a restricted A^-single- 

body method introduced in |2Zj. Here, test particles sampled from the initial HDM phase 

space distribution are allowed to evolve, one at a time, in an external potential </)(x, r) 

due to the CDM halo alone. Gravitational interactions between the HDM particles 

themselves are explicitly ignored, which is well justified considering current cosmological 

constraints on the HDM mass fraction today (^HDM/^m ;$ 0.1). Compared with other 

solution methods, our A^-single-body approach is computationally inexpensive relative 

to full-scale A^-body simulations, and is yet able to reproduce essential nonlinear effects 

not captured by simpler, linear methods (e.g., [52!). See reference [2Zj for more details. 



Appendix B. The halo model 

The halo model supposes that all matter in the universe is partitioned into distinct 
units, the halos. This assumption allows one to study the universal matter distribution 
in two steps: the distribution of matter within each halo, and the spatial distribution 
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of the halos themselves. In its simplest application, the halo model proposes that the 
matter power spectrum P(k, z) be composed of two distinct terms (e.g., [36,), 

P(k, z) = P^^ik, z) + P 2 - halo (fc, z), (B.l) 



where 



P 1 -^ 10 ^, z) = J dM n(M, z) 



' M 

, Phaio, 



> 2 - halo {k,z) = P lin (k,z) 




dM n(M,z) b(M) - — ] u(k\M) 

V Phaio 



(B.2) 



describe, respectively, the correlations between two points drawn from the same halo 
("1-halo") and from two different halos ("2-halo"). The former dominates on small 
scales (i.e., large fc's), while the latter rises to prominence on large scales (i.e., small fc's) 
and approaches the power spectrum calculated from linear perturbation theory P hn (k, z) 
as k — > 0. The average (comoving) matter density phaio counts all matter clustered in 
halos. For a basic ACDM set-up, phaio is well approximated by phaio — Pm = ^mPcrit, 
where p crit is the present critical density. 

Three pieces of information are required to complete the model. 

(i) The mass function n(M, z) specifies the comoving number density of halos of virial 
mass M at redshift z. Following from the Press-Schechter theory |HS1, this is usually 
written as 

n (M,z)dM=^f(u)du, (B.3) 

where f(y) is a universal function of the peak height v = S^ c /af in (M, z). Here, 
5 SC = 1.68 is the linear overdensity at the epoch of spherical collapse, and 

o-L(Afz) = J ^P^(k,z)\W(kR)\ 2 (B.4) 

is the rms of the linear fluctuations filtered with a tophat window function 
W(x) = 3/x 3 (sinx — xcosx) on a scale of R= (3M/47TP™) 1 / 3 . In this way, every 
halo mass M is associated with a unique value of u, with large i/'s corresponding 
large halo masses and so forth. A number of variants for the function f(y) exists 
in the literature. Here, we adopt the version proposed by Sheth and Tormen 



vf{v) = A (l + (g„)-*) Jf.exp{-qu/2), (B.5) 
v ' V 271 

with p = 0.3 and q = 0.707, based on fits to A^-body simulations. The constant A 
is determined by mass conservation, i.e., / dM n(M, z) M = phaio, or, equivalently, 

/£*!//(!/) = 1. ^ 

The linear bias b(y) parameterises the clustering strength of halos relative to the 
underlying dark matter, and obeys / dv f[y) b{y) = 1 by construction. For the 
Sheth-Tormen mass function (jB.5|) . the associated bias is 



K„) = 1 + ^ + (B.6) 
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(iii) The function u(k\M) is the Fourier transform at wavenumber k of the matter 
distribution p(r\M) within a halo of mass M, normalised to the halo mass. For 
spherically symmetric density profiles p{r\M) truncated at the virial radius r vir (to 
be defined below), this is given by 

,,,„, p(k\M) rvir 2 sin At p(r\M) 

u(k M) = HK i- ' = / dr 4vrr 2 — Hy ' ' . B.7 

M Jo kr M v ; 

In the simplest case where only CDM and baryons cluster, a natural choice for 

p(r\M) is the Navarro-Frenk- White (NFW) profile 

p(r\M) = Ps , (B.8) 

{r/r s )(l + r/r s ) 2 

which has proven to provide a very good description of the density run around 
virialised halos in high resolution iV-body simulations in a variety of CDM 
cosmologies. Here, r is the radial distance from the halo centre. The parameters r s 
and p s are determined by the halo's virial mass M and a dimensionless concentration 
c = r vir /r s via 

_ ^th cr* 

Ps 3 ln(l + c)-c/(l + c)' 1 ' 

r * l l 3 M V /3 minl 

r s = = - -T^— , ( B - 10 

C C y 4:710th Pm ) 

where 5th is the virial overdensity predicted by the spherical top-hat collapse model 



TH ^ 



18tt 2 + 82y - 39y 2 



"""W" 1 ' n w = d^- (R11) 

In addition, a tight correlation between c and M has been observed in the analysis 

of [SE], 

9 /M\~ - 13 

^ ~ T+~z (if) • < R12 » 

where M* = M{y = 1) denotes the P lm -dependent nonlinear mass. 

Several modifications to the halo model have been discussed in the literature, to 
account for, e.g., the effects of baryons neutrino infall [37], etc.. To include the 

effects of HDM infall, we follow the prescription of [H7j by defining 

M = M CDM+ 6 + Mhdm, (B.13) 

u(k\M) ee /Wfc|M) + p CDM+fe (/cW (R14) 
McDM+b + Mhdm 

where Mj and pi{k\M) denote, respectively, contributions from the i type particles to the 
total halo mass and Fourier transform of the density profile. Furthermore, the quantity 
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Phaio should account also for HDM residing in halos, i.e., phaio = Pcdm + Pb + PHDM,haio — 

(^CDM + ^fe)Pcrit + PHDM,halo, where 

<OHDM,halo = Phalo J dvf{v) ™ M ■ (B.15) 

Note that, in general, pHDM.haio "C ^HDMPcrit, since the HDM's significant thermal energy- 
prevents a large fraction of particles from participating in the infall. 
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